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Abstract 

The last stable orbit (LSO) of a compact object (CO) is an important 
boundary condition when performing numerical analysis of orbit evolu- 
tion. Although the LSO is already well understood for the case where 
a test-particle is in an elliptical orbit around a Schwarzschild black hole 
(SBH) and for the case of a circular orbit about a Kerr black hole (KBH) 
of normalised spin, S (\J\/M , where J is the spin angular momentum 
of the KBH); it is worthwhile to extend our knowledge to include ellipti- 
cal orbits about a KBH. This extension helps to lay the foundation for a 
better understanding of gravitational wave (GW) emission. 

The mathematical developments described in this work sprang from 
the use of an effective potential (V) derived from the Kerr metric, which 
encapsulates the Lense-Thirring precession. That allowed us to develop a 
new form of analytical expression to calculate the LSO Radius for circular 
orbits (Rlso) of arbitrary KBH spin. We were then able to construct a 
numerical method to calculate the latus rectum (Ilso) for an elliptical 
LSO. 

Formulae for E 2 (square of normalised orbital energy) and L 2 (square 
of normalised orbital angular momentum) in terms of eccentricity, e, and 
latus rectum, I, were previously developed by others for elliptical orbits 
around an SBH and then extended to the KBH case; we used these results 
to generalise our analytical Ilso equations to elliptical orbits. LSO data 
calculated from our analytical equations and numerical procedures, and 
those previously published, are then compared and found to be in excellent 
agreement. 

04.20.-q, 95.30.Sf, 95.75.Pq 
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1 Introduction 



One of the most important goals in experimental gravitation today is the de- 
tection of gravitational wave (GW) radiation [TJ |2] EJJ. To achieve this goal, 
considerable effort has been made to improve the theoretical understanding of 
the evolution of compact object (CO) orbits in extreme black hole systems 
jU [5] |6l El HI |9] • In this paper, we assume a point-like test-particle since the 
ratio of CO mass to the mass of the massive black hole (MBH) will be small 
(i.e. < 10 HO]); and the effect of the CO mass upon the post Newtonian 
(PN) equations that we will use in our subsequent modelling of the CO orbits is 
negligible |llll2lfl]. In this paper, we shall then assume the behaviour of the CO 
to be closely approximated by that of a test-particle. If reference is made to the 
orbital evolution of a true CO, as described by the PN evolution equations, then 
we will use the term, CO. The treatment of CO orbital evolution we will present 
in a forthcoming paper will be based on the work of [4] |5| [T5J HU E] 13 [HI US] in 
which PN equations for a rotating MBH, also called a Kerr black hole (KBH), 
are considered. 

The objective of this study is to lay the foundation for our subsequent work 
that will include the numerical calculation of the GW energy emission by ex- 
treme KBH systems where the CO is in an elliptical orbit in the equatorial plane 
of the KBH. The most basic quadrupole model |4j [5] admitted solutions in closed 
form [5][TJ]; but because the more comprehensive evolution equations now used 
are too complicated to admit an analytical solution, numerical integration of or- 
bital parameters is required [8]. Therefore the last stable orbit (LSO) becomes 
an important boundary condition. Such an analysis requires an understanding 
of how V depends on L and on the inclination of the CO orbit. To undertake 
future work for inclined orbits it is important to know the minimum physically 
meaningful value of L. 

Previous research has demonstrated how the effective potential (V") of a 
test-particle in an equatorial orbit around a Schwarzschild black hole (SBH) 
|1 51 [T6j IT?] can be calculated from the Schwarzschild metric and used to deter- 
mine the latus rectum of the LSO (Ilso)- A treatment of V for a KBH system, 
where the test-particle follows a circular LSO (section 12.7 of [18]) yields an an- 
alytical expression for the value Rlso m terms of normalised spin, S, (equation 
12.7.24 in [18], [19]) (S = \s\/M where s = 3/M and J represents the spin angu- 
lar momentum of the KBH). Such treatment of V also gives rise to expressions 
(equations 12.7.17 and 12.7.18 in [18]) for the orbital energy, E, and orbital 
angular momentum, L. In |16] the energy and orbital angular momentum equa- 
tions were also derived for an SBH system with the test-particle in an elliptical 
orbit. In the significant work by Glampedakis and Kennefick [2], their treat- 
ment of E and the quantity, (^L — SEJ , enabled us to derive generalised Rlso 
formulae for elliptical orbits. Analytical expressions have a clear usefulness in 
the development of new theoretical concepts and numerical methods |10l [20] . 

The Lense-Thirring effect, an apodeictic pTj prediction of general relativity, 
is the means by which the rotation of the KBH imparts important changes on 
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the test-particle orbit |22l |23l |2"T1 [23] that are distinct from those associated 
with the SBH. The swirling of spacetime in the vicinity of the KBH applies a 
torsion to the orbiting test-particle; therefore, the orbit evolution will be altered, 
thus causing changes in the point at which the test-particle reaches its LSO. We 
shall develop an analytical and numerical methodology to calculate the LSO of 
a test-particle in elliptical orbit about a KBH. Numerical estimates of the latus 
rectum of the elliptical LSO orbits with respect to KBH spin are available in 
the literature (Table I in (2], based upon the work of Schmidt [25], and Table I 
in [10]); and they will provide a means to validate our results. 

In section [27T] the Kerr metric is introduced and used in section [272] as the 
basis of developing some essential analytical formulae to calculate the orbital 
angular momentum of test-particles in circular paths around a KBH (section 



2.3 1. A formula for Rlso (prograde and retrograde) is then developed analyti- 



cally; and the general formula for the Ilso of elliptical orbits is also presented. 
In section [274] the development and demonstration of a numerical algorithm to 
determine the latus rectum and eccentricity of test-particles of higher orbital 
angular momentum then follows. The results of this analysis (section |2.5| , as 
well as results obtained from the general analytical formulae for LSO latus rec- 
tum, are compared with results obtained from the literature. In section [3] we 
shall draw conclusions. 



2 Understanding the Last Stable Orbit About a 
Rotating Massive Black Hole 

2.1 Kerr Metric 

The Kerr metric (See equation 13.12 in [26 ) represents the solution to the 
Einstein Field Equations in the case where the MBH possesses spin angular 
momentum, 
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in which the 



factors, R = r/M and S = \3\/M 2 , are used to express the metric in di- 
mensionless terms. The symmetric off-diagonal elements, — 2M 3 SR sin 2 (9) / p 2 , 
correspond to the Lense-Thirring precession that arises from the spin of a cen- 
tral KBH of mass, M. Observe that when, S = 0, the Kerr metric equals the 
Schwarzschild Metric. 
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Although the Schwarzschild Metric is expressed in spherical coordinates, 
when the central black hole rotates it is appropriate to use the Kerr metric 
expressed in Boyer-Lindquist (BL) coordinates. The conversion of the BL co- 
ordinate system variables to Cartesian coordinate variables is represented by 
these equations (see equation 11.4.7 in [27] and see also [28]) : 



x = \J r2 bl + S 2 sin (<?) cos (0 - /) j 

V = \/R% L +&sin(e) S m(ct>-f), 

z = R BL cos(9), (2) 

where 

/ = iarctan {s/r B l)- (3) 

Because < S < 1.0, a prograde or retrograde orbit is represented by the 
respective use of a plus or minus sign in equation The BL coordinate 

system will be used throughout this treatment. The conversion of LSO radius 
from BL to spherical coordinates is required whenever one performs a simulation 
of the evolution equations reported in |29l [7J [30] |3T] |8]. This conversion is 
uncomplicated in the current application (in which the angle, t, between the 
orbital angular momentum vector and the spin axis of the KBH is zero), and 
proceeds by adding the squares of x, y, and z as shown in equation (|2j to obtain, 

^Spherical — X 2 + y 2 + Z 2 . (4) 

By substituting the relationships in equation ^ into equation Q , one obtains 
the mathematical relationship, 

^Spherical = R BL + S^SHl 2 (9). (5) 

Recall that S is the normalised spin of the KBH and 8 is the polar angle of the 
test-particle in its orbit. In this study, we work with orbits that are exclusively 
in the equatorial plane of the KBH. Therefore one sets 6 — ^ to obtain 

^Spherical = R BL + (6) 

Such a relationship is required for transforming LSO radii (BL coordinates) into 
the spherical coordinate system. 



2.2 Effective Potential 



We shall develop a formulation of the effective potential of a test-particle in 
orbit about a KBH. By so doing, the location of the LSO can be estimated. 
In the following equations and calculations the radius, R, is represented in BL 
coordinates. For simplicity of notation, the BL subscript will be suppressed 
(except in Section 2.5.21. 
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The four-momentum can be expressed as: 

- 2 'dR 



P, = 
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for zero mass, where E is the energy of the orbital element and L is the or- 
bital angular momentum of the particle in orbit normalised with respect to its 
mass, to, and the KBH mass, M. The (dR/dr) is the derivative of the radial 
component of the compact object with respect to the proper time, r. For the 
zero-mass particle (which has no rest mass), L is its total linear momentum (vis. 
L = E. p i lo ton/c). The factor, dR/dX, is the derivative of the radial component 
of the zero-mass particle with respect to an affme parameter, A, which is used 
in place of proper time, r, since a zero-mass particle always follows a null path. 

The invariant quantity of mass-energy can be calculated for each case of a 
test-particle of infinitesimal mass 



P P 

and a zero-mass particle 



#7 I 



I Kerr 



(9) 



PP = P y P sS ^\ Karr = 0. (10) 

In that respect, the expected behaviour of a test mass will differ from that 
of a zero-mass orbital element. From these equations, the effective potential 
can be calculated by making a few assumptions about the path taken by the 
orbiting zero-mass particle. The inverse Kerr metric (g 61 ) is shown in Appendix 



A (equations (77 1 and (781) 



2.2.1 Test Particle. 

We restrict our work to the case of a test particle of mass, to, in orbit about a 
KBH with 9 — |. By evaluating P ■ P (see equation (j9j) using the test mass 
four-momentum (see equation (|7]l) one obtains, 

P-P =- (^R 4 E 2 -R 4 m 2 (^j 2 ~R 2 m 2 L 2 + R 2 E 2 S 2 
+2 RE 2 S 2 + 2 Rm 2 L 2 + 4 RESmlj 



x (i? 4 -2i? 3 + i? 2 5 2 ) 



-l 



,2 



(11) 
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To develop a relationship between the effective potential and the orbital 
para meters several sequential steps must be followed. First, all terms in equation 



(11 1 are collected and equated to zero, and then divided by m 2 and the (dR/dr) 2 
terms are collected on the right hand side of the equation. Noting that E/m — E 
represents the specific energy content of the orbiting test-particle, one then 
obtains, 



[r 2 S 2 + 2RS 2 + i? 4 ) E 2 - (aRSL) E 
[L 2 R 2 - 2L 2 R + R 2 (R 2 -2R+ S 2 ^ ) 



dR" 2 



At the points of closest (pericentre) and farthest (apocentre) approach the 
derivative of R with respect to t is zero. By performing that simplification, 
one obtains a quadratic equation in E, i.e. 

- ( R 2 S 2 + 2RS 2 + i? 4 ) E 2 + UrSL\ E 

L 2 R 2 - 2L 2 R + R 2 (r 2 -2R+ S 2 ^ 
= 0. (13) 



The factored form of equation (13 1 corresponds to the following equation |17] : 

(E - V+^j (E - V-^j = 0. (14) 
Therefore two solutions for the effective potential can be calculated: 



a 
b 



—b =p Vb 2 — Aac 
2a 

- (i? 4 + R 2 S 2 + 2 S 2 R^j 
ARSL 

L 2 R 2 - 2L 2 R +R 2 (r 2 -2R + S 2 



(15) 



For the SBH (i.e. S — 0), the value of V± (from equation (15)) becomes: 



V 2 = 



(R - 2) (r 2 + I 2 ) 



i? 3 



(16) 
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which depends only on the values of R and L, as expected (as shown in Figure 
0- 

The effective potential contains important information. In the case of the 
SBH, the relationship between V± and R describes the test-particle orbit and 
leads us to a calculation of the values of L and R at which the test-particle can 
no longer sustain a stable orbit. The LSO is an important characteristic of the 
binary system that is identified as the point at which the V + curve (Figure]!} has 
a slope of zero and the second derivative with respect to R is not positive. The 
effective potential, V-, corresponds to particles and photons for which their 
orbital angular momentum has an opposite sense to the KBH spin (section 
11.3 in |17| ). The mathematical treatment of V+ presented in the sections 
that follow preserves its prograde and retrograde properties; indeed, we have 
found that the use of V- in the calculations that follow yield the same results. 

2.3 Last Stable Orbit (LSO) for a CO in the Equatorial 
Plane of the Kerr Black Hole 

The equations for the radius of a circular or elliptical LSO can be calculated 
through a mathematical treatment of the following two equations: 

§ - » 

and 



£ <- »■ 

where the point of inflection (which corresponds to a circular LSO) can be found 



by evaluating the intersection points of the equations (17 1 and (18 1. 

The loci of these two equations is depicted in the \~R, plane for a KBH 

with a spin value of S = 0.5 (see Figure Their intersection points (de- 
rived numerically with Maple 11), [ R = 7.554584715, L = -3.884212633 ] 
and [ R = 4.233002530, L = 2.902866150 ] , correspond to the radial position 
of the LSO, R, of a test-particle with an orbital angular momentum of L. These 
points differ from [ R = 6.0, L = ±V12 1, which is the solution for an SBH. 
The existence of an intersection point on the graphical plot notwithstanding 
(see Figure [2}, on frequent occasions, no result was returned by Maple. On 
other occasions a correct value of R was returned, while the value calculated for 
L deviated by at least a factor of two from the graphical result. Such inconsis- 
tent behaviour was attributed to the great complexity of the expressions being 
treated and the associated floating point round off error; therefore, an analytic 
method was sought. 
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2.3.1 Orbital Angular Momentum. 

The derivative of V+, equated to zero, can be used to determine an analytical 
expression for L 2 in terms of R and S for circular or elliptical orbits. From 



equations (15 1 and (17) one obtains 



dV+ 
~dR 



= \^-3S 2 R 4 L 2 + 6S 2 R 3 L 2 -2R 2 S 4 L 2 

+R 5 L 2 S 2 + 2 S 2 R 5 -3R 6 S 2 - 3 R 6 L 2 

-3 S 4 R 4 - S 6 R 2 -2RS 6 + 8 S 4 R 2 - R s + R 7 L 2 



+ (6R 2 SL + 2 S 3 L) Vi? 3 (R 2 - 2 R + S 2 ) (i? 3 + L 2 R +S 2 R+2 S 2 )^j 

+ (^/R 3 (R 2 -2R+ S 2 ) (i? 3 + L 2 R + S 2 R + 2 S 2 ) 

(R 3 + S 2 R+2S 2 f^j 
= 0. (19) 



The denominator of equation (19 1 can be disregarded because the quotient is 
equated to zero; it is also required that the roots of the factors present in 
the denominator lie outside the range of physically attainable R values. To 
be specific, the roots of ^R 2 — 2R + S 2 ^j correspond to the event horizon for 
massless particles, those of R 3 are zero and beyond the LSO, and the roots of 
(i? 3 + L 2 R + S 2 R + 2 S 2 ) and (i? 3 + S 2 R + 2 S 2 ) are complex and thus also 
physically unattainable for real values of R. 

The simplified power series is thus derived from the numerator of equation 



(19) after eliminating the square root, 



-2R 2 [ -i 



R 3 [9R-6R 2 + R 3 ~AS 2 jL 4 

(-3 R 4 + R 5 - 12 S 2 R + 6 R 2 S 2 + 2 R 3 S 2 + 5 S 4 + S 4 R^ L 2 

'^R 4 + 2 R 2 S 2 — 4 S 2 R + S 4 ^j 2 
= 0. (20) 
Therefore L 2 can be obtained directly by using the quadratic formula, 



-b ± \/b 2 - Aac 



L Z = ~ , (21) 

la 



where we have redefined: 



R 3 (9 R - 6 R 2 + R 3 - 4 S 2 ) 
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b = 2 If [-3 i? 4 + R> 12 S 2 R + 6RW + 2 + 5S* + S*R) 

i 



c = R 



2R Z S 2 -4S Z R+ S' 



Two solutions are found that correspond to the orbital angular momenta of a 
test-particle in a prograde orbit, 
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-3 R 6 + R 7 — 12 R 3 S 2 + 6 R 4 S 2 

2R 5 S 2 + 5S* 4 i? 2 + S 4 R 3 

2 S (3 R 2 + S 2 ) (r 2 -2R + S 2 ) VR 3 ^ 



(r 3 (9R-6R 2 + R 3 -AS 2 
and in a retrograde orbit, 



(22) 



J Ret 



3 i? 6 + R 7 - 12 R 3 S 2 + 6 i? 4 ^ 2 

+2i? 5 5 2 + 5 5 4 i? 2 + S' 4 i? 3 
+2 £ (3 i? 2 + 5 2 ) (i? 2 - 2 + S 2 ) \ /? ' 

i? 3 f 9 i? - 6 R 2 + R 3 - 4 S 2 



(23) 



An analytical expression for L 2 with respect to R and 5 has been derived. 
But one must consider that the formula is limited to providing a value of L 2 
that corresponds to a test-particle in its LSQ (BL coordinates) about a KBH of 
spin S. These formulae (equations (22 1 and (23 1) do not provide a relationship 



between R and L 2 for a general orbit. 

Consider an example where S = 0.5. The relationship between the value of L 
and the radius R € [1.0, 6.0] is plotted in Figure|3] One observes a power series 
for which the values of R and L for a circular orbit (at the point of inflection) 
occur at the local minimum. Therefore it is possible to derive an expression for 
the radius of the LSO, Rlso, a t that point of inflection for an arbitrary spin, 
S, where < S < 1. 

2.3.2 Circular LSO Radius. 

The calculation of such an analytical relationship proceeds as follows. The 
derivative of L 2 with respect to R is set equal to zero. Prom equation (21 1 we 
obtain: 



d(l 2 
dR 



SVR/R ( 3R 7 - 45i? 5 + 20S' 2 i? 5 + 54i? 4 - 265 2 i? 4 
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-9S 4 i? 3 + 24S 2 i? 3 - 26S 4 R 2 - 54S 2 R 2 + 53S 4 R - 12S 6 ) 



± 



(i? 8 - 2i? 6 S 2 - 3R 4 S 4 - 12R 7 - 28S 2 R 5 - 24S 4 R 3 
+45i? 6 + 126S 2 i? 4 + 57S 4 i? 2 + 20S 6 

-54i? 5 - U4S 2 R 3 - 90S 4 R+ 108S 2 R 2 



R 2 



-R 3 + AS 2 + 6R 2 - 9R 



0. 



(24) 



Where the plus sign corresponds to a prograde orbit and the minus sign corre- 
sponds to a retrograde orbit. The denominator contains a factor (i.e. R 3 + AS 2 

with roots that correspond to the event horizon and a factor R 2 with roots equal 
to zero, which lie beyond the event horizon and are thus unattainable. 

The simplification of the equation by taking only the numerator and elimi- 
nating the square root, can proceed to yield the following result: 



6R 2 



R 3 ( S 2 - R J [9 S* 4 - 28 S Z R - 6 S 2 R 2 + 36 R 2 - 12 R 3 + R 

x (S 4 + 25 2 i? 2 - 4 S* 2 i? + i? 4 ) (r 3 -4S 2 -6R 2 + 9R?) 2 

= 0. (25) 

Fortunately, the polynomial that expresses the relationship between S and 
Rlso, is already simplified into a product of some binomials, trinomials, and 
quartics (see Table [TJ. Each one can be assessed by considering the examples of 
an SBH with no spin (S — 0.0) and a KBH with S — 0.5. For the former case, 
the solution, [ R = 6.0, L = \fV2 ] , is known; for the second case, it has been 
calculated numerically, [ R = 4.233002530, 1 = 2.902866150 ]. These cases 
help one to identify the relevant factor. It is interesting to observe that some of 
the radii in Table [l] have complex values. 

The factor that yields the values of the LSO radii (one for each of the possible 
prograde and retrograde orbits of the CO) is: 



9 S 4 - 28 S 2 R - 6 S 2 R 2 + 36 R 2 - 12 R 3 + R 4 = 



(26) 



This quartic equation ( 26 1 can be converted to a companion matrix which is 



solved for its eigenvalues to yield the analytical solutions for Rlso for the 
prograde and retrograde orbits (see Appendix B). These solutions are: 
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and 




Z + 3 (3 + S 2 ) (28) 



where: 



Z = 3 + S 2 + (3 + S) + 

+(3-5) ((1-S)(1 + S) 2 )\ 



Although formulae that are analytically the same as ours have already been 
developed by Bardeen et al. [19], our formulae were derived by independent 
means and are simpler. The numerical results of each equation differ insignif- 
icantly over the physically valid range of < S < 1.0. And our formulae are 
more robust with respect to round-off error when evaluated numerically, and 
they are roborant of the pre-existing calculations. 

2.3.3 Orbital Energy and Angular Momentum at the LSO. 

One can derive new formulae for the test-particle orbital energy, E, and angular 
momentum, L, in terms of parameters S, e, and latus rectum, I, by using 



equation (15 I as a starting point. We know that, 

V+ = E, (29) 
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2RSL- 



R R 2 - 2R+ S 2 ) R 5 + R 3 L 2 + R 3 S 2 + 2R 2 S 2 



R[R 3 + S 2 R + 25' 



(30) 



when the test-particle is in its LSO. Although the roots in R are readily found 
by Maple, they are inordinately long and not useful. A more effective derivation 
method, similar to the one used in [2], shall be outlined. 
By manipulating the formula in equation (30 1 we obtain: 



R s - 
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(31) 
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from which we can obtain the expressions for the sum and the product of the 
roots in R directly from the coefficients of the polynomial, vis. 

(R- ri )(R-r 2 )(R-r 3 ) = (32) 

which implies, 

R 3 -(r 1 +r 2 + r 3 )R 2 + (r 1 r 3 + r 1 r 2 + r 2 r 3 )R-r 1 r 2 r 3 = 0. (33) 



Where {n , r 2 , r 3 } are the roots in (|32| and ( 33 1 . We find the following equations 



for the sum of the R roots (i.e. R sum — n + r 2 + r 3 ): 

Rsum = 2 (l - E 2 ^' 1 , (34) 
and for their product (i.e. R pro d = fir 2 r 3 ), 

Rprod — Rsum {l — ES\ . (35) 

The corresponding formulae for E and L are as follows: 

_j_ V Rsum (Rsum ~ 2) (36) 

Rsum 

and, 

j \J Rsum (Rsum 2)tS i ^/ ' RsumRprod /otA 

For the LSO, the roots, {ri,r 3 }, correspond to the LSO radius; therefore, we 
make the following substitutions: 



ri = r 3 = R Mm = — — , (38) 
1 + e 

and 

I . . 

r 2 = Rmox = ; , (39) 

1 — £ 

where / is the latus rectum of the elliptical LSO. We can now set: 



Rsum — %RMin ~t"" Rlvlax 
_ I I 



l + e-l-e (40) 



and 



13 



Rprod — RMin R-Max 
13 



(l + £ y(l-e) 



(41) 



By substituting equations (40 1 and Mil into equations (36 1 and (37 1 the 



following formulae are obtained: 



and, 



E = 1-2 



(42) 



L = \ SE ± h 



(1 + e) (3-e) 



(43) 



fL-S£j) (1 + e) (3-e) 



= i 2 



(44) 



They express the square of the orbital energy and the orbital angular momentum 
in terms of the eccentricity, e, and latus rectum, I, of a test-particle in its LSO 
about a KBH of spin, S. In equation (43 1, the prograde orbit takes the minus 



sign and the retrograde orbit takes the plus sign. The modified form of equation 



(43) shown in equation (44 1 corresponds to equation (23) in j2]. 



Similar equations derived by Cutler, Kennefick, and Poisson [16] 



E 2 = 



((J- 2 (1 + £)) (J- 2 (1 - e))^ r 1 (J -3- t 



(45) 



and 



I 2 1-3 



(46) 



are only valid for SBH systems. Equation (43 1 reduces to equation (46 1 when 
S = and the relationship I = 6 + 2e is used. 

Glampedakis and Kennefick [2] present a similar treatment which has the 
advantage of yielding more general results since it is not assumed that the test- 
particle has reached the LSO (i.e. r\ > r^j. Therefore 



r-3 



2(Z-^) 2 (l-e 2 ) [F(l-£ 2 ) 



(47) 
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with, n = i?Mm and T2 = Rmo,x, as before. Their formula for energy, 
E = 



'l - F 1 (i - e 2 ) ji - F 2 (Z - (1 



(48) 



proves to be ideal for generalising our formulae for circular LSOs, Rlso, to one 
for elliptical orbits, Ilso (See Appendix C). 

2.3.4 Elliptical LSO Radius. 

The evaluation of X 2 = (^L — SEJ in |2 provides a means to extend equations 
(27) and (28) beyond their use with circular LSOs to more general elliptical 



LSOs by direct substitution of X into equation (44 1. Although a leading order 



Taylor expansion (see equation (24) in [2] ) is available from a slow rotation 
approximation of equation (44 1 (i.e. S ~ 0); we present our analytical results. 



The general form of the Ilso equations for elliptical orbits are: 



and 



Where: 



(3 + e) + y% 



" 1.6 S2 ^J. £) -Z + (3 + ef + S 2 (l + s) (3 - e) 



(3 + s) + V^o 
S 2 (l+s) 



16 



Z + (3 + ef + S 2 (1 + e) (3 - e) 



(49) 



(50) 



l/3S 2 (l + e) (3 - e) + 1/3 (3 + e) 2 

S 4 (1 + ef (3 - ef - 2 S 2 (3 + e) (1 + e) (e 2 + 15) + (3 + ef 



1/3 



1/3 (^) 3 , 



(Zi) 



Zi = (3 + e) 6 

+ S 2 {l + e) (5 2 (l + e) (s ,2 (l + e)(3-e) 3 + 3e 4 + 18e 2 +459) - 3 (e 2 + 15) (3 + e) 3 ) 
+ 24 VS^Zu, 

and 

Z« = (l + e) 4 ^ 6 (l-^ 2 ) ((l-e)(e + 3) 3 -5 2 (l + e)(3-e) 3 ) . 
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As required, equations (491 and (50 1 reduce to equations (27 1 and (281 when 
£ = 0. By setting S = 0, both equations reduce to I — 6 + 2e. And in the 
extreme cases, where S — 1 (retrograde and prograde), equation (49 1 reduces to 
(]50| reduces to 5 + e + 4\/l + e , as required. A detailed 



I = 1 + s; and equation 
treatment of equations (|49 



and (501 will be outlined in a forthcoming paper. 



2.4 Calculating the LSO Properties 

2.4.1 Introduction. 

The elements have now been found to perform general calculations of the LSO 
for arbitrary values of KBH spin, S, orbital angular momentum, L, and total 
energy, E. Although we have analytical formulae that give us Ilso f° r general 
elliptical orbits, it is important to construct and outline our methodology in 
preparation for future work on test-particle orbits that are inclined with respect 
to the equatorial plane of the KBH. We must quantify the relationship between 
the value of L for the test-particle orbit and the shape of its effective potential 
surface. 

Here we outline, in detail, our numerical method for calculating the latus 
rectum, Z, and eccentricity, e, of LSO orbits. These values will help us to appraise 



the usefulness of our new, generalised Ilso equations in (49 1 and (50 1 



2.4.2 Algorithm. 

For clarity, an example where S = 0.5 and the test-particle is in a prograde 
orbit is demonstrated. In Table [2] the calculations for a retrograde orbit, and 
an SBH are included for comparison. 

Such an algorithm proceeds as follows: 



Specify the KBH spin - A given problem will most likely have a prior 
specification of a fixed value of S, where < S < 1 for either a prograde or 
retrograde orbit (if a retrograde orbit is used, (ret), will follow the value assigned 
to S). In this example we shall use S — 0.5 and a prograde orbit since it has 
already been used in the calculation of Rlso an d L for prograde and retrograde 
LSOs previously in this paper (see 
performed for the S = 0.5(ret) case, and for an SBH (See Table [2j. 

We use either equation ( |27| for a prograde orbit or equation ([28} for a 
retrograde orbit to directly calculate Rlso (BL coordinates) thus, 

S = 0.5 => Rlso = 4.23300, (51) 
which gives us the LSO radius of a circular orbit, Rlso- 

Find L - The values of S and Rlso can now be used to calculate the value 
of L 2 assuming the LSO is at a point of inflection (i.e. a circular orbit) vis. 



2.3 and Figure^}. Similar calculations were 
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equations (22 1 or (23 1 depending on the direction of the orbit. The result for 
S = 0.5 and R LSO = 4.23300 is found to be, 



L 2 = 8.4266319. (52) 

This value is necessarily a positive quantity, hence the need to ensure that the 
correct prograde or retrograde orbital angular momentum equation has been 
used. 

Calculate E - Because the values of S, Rr,p n, and L are known at the point 



of inflection, we can use V+ (see equation (15 l) to directly calculate the energy, 



E, of the test-particle in a circular orbit, i.e. 



E = 0.91788201. (53) 

N.B., the value of E < 1.0, hence the orbit is bound. Whenever E ^ 1.0, the 
orbit is not bound. 

Expand to include elliptical orbits - By careful examination of Figure [3] 
one sees that the local minimum of L corresponds to the case where the LSO 
is circular; the values of Rlso an d the radius of the local minimum of the 
potential, V+, (Figure [lj coincide, as expected. The angular momentum, L, 
that corresponds to an elliptical LSO is then higher than that for a circular 
LSO. 

The algorithm shall be broadened to include the case of an elliptical orbit. 
For the orbit to be elliptical, the orbital angular momentum, (L — yip- => L = 
V8. 4266319) must be increased by an arbitrary factor SL (where SL > 0, see 
Figure [i}; the slight increase in L above its minimum value changes the LSO 
from a circular orbit, to one that is elliptical. Accordingly the value of E will 
increase and the value of Ruin will be reduced. A similar treatment of elliptical 
orbits, based upon increments of L, can be found in jlOJ . 

Find Rmui for the elliptical orbit - By working with a larger value of 
orbital angular momentum in the form (L Elliptical — ^Circular + SL) we can 
calculate the value of Rutin without requiring the new value of the orbital energy, 
E (see Figure If SL = 0.01, then L = 2.904588078; therefore, (vis. equations 



(22| or (p3f) the new value of Ruin can be calculated numerically to yield: 



R Min = 4.10329200. (54) 
Correspondingly, the total orbital energy can be calculated (vis. equation (|T5|): 



^Emptied = 09180746j (55) 
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^Circular = Q.91788201. (56) 

As required: Ef^ Ucal > E^ ular . 

Find the maximum radius for an elliptical orbit - In calculating a data 
set, the various values of 5L are selected and the corresponding values of e and 
fare found. Now that the value of ^ l J^P tlcal [ s known, the maximum radius of 
the elliptical orbit (R,Max) can be calculated numerically, vis. V+ = E LS Q Uca , 
because the effective potential of the test-particle has the same value at Ruin 
and Rmo,x- The result is: 

Rmcx = 4.520999771. (57) 

Determine the elliptical orbit parameters - We now have the information 
necessary to calculate the latus rectum, I, and the eccentricity, e, of the orbit. 
The dimensionless semi-major axis, A, of the elliptical orbit can be calculated 
from the values of Rlso an d Rmux'- 

. (RMin + RMax) 

A ~ 2 ' 

= 4.3121. (58) 



= 0.0484. (59) 
The latus rectum, /, is now be calculated, 



I = A(l-e 2 ), 

= 4.3020. (60) 

N.B., the values of A, I, Rlso (f° r circular LSO), Ruin and Ruax are expressed 
in terms of BL coordinates. 

2.4.3 Calculation of the normalised orbital frequency 

According to the relativistic form of Kepler's third law (see problem 17.4 in 
|32] or exercise 12.7 in [18 ), the orbital period of a closed orbit, P, can be 
expressed in terms of the semi-major axis of the orbit, a, and the mass, M of 
the central body about which the test-particle orbits. This equation applies to 
elliptical orbits in general. If the orbit is subject to precession, then the value, 
v, represents the orbital frequency of the test-particle in an open orbit. Hence, 
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where s = 3/M; and the plus sign corresponds to the prograde orbit and the 
minus sign corresponds to the retrograde orbit . The corresponding orbital 
frequency is: 



v = P~ 



(62) 



therefore, 



1 

2^ 



M 



(a 3 / 2 ± |s| y/Ht) 



(63) 



which leads to, 



W = 2nMv 



M 3 / 2 



(a 3 / 2 ± |s| 



'M 



(64) 



(in equation (64 1 the parameter, o, refers to the length of the semi-major axis). 
When variables normalised with respect to the KBH mass, M, are used: 



and 



one obtains, 



a 

M ' 



|s| /M; 



(65) 



(66) 



W = 



(67) 



If equation (60) is then used to represent equation (67 1 in terms of the dimen 



sionless latus rectum, I: 



W 



(l-£ 2 )7(Zl±S(l-£ 



2\ 2 



(68) 
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2.5 Calculations 



2.5.1 LSO and orbit characteristics. 



Three methods were used to calculate Ilso f° r orbits of various eccentricity 
(0 < £ < 1.0) and KBH spin (S = 0.5,0.99; prograde and retrograde). The 
values obtained here are shown alongside those found in the literature |10l [2] 
m Tables ^ and ^ The Ilso values we obtained by following the algorithm 
described in |2.4.2| are listed in the Numerical column. The general formulae 



described in |2.3.4 were used to generate the values in the Analytical column. 



A third method was used to numerically estimate the Ilso values directly from 
the companion matrix (Appendix C) by first substituting the S and e values 
into the matrix before calculating its eigenvalues. 

The agreement between our various calculation methods, and with the re- 
sults published previously in |101 12] is excellent (i.e. error < 0.1%). Therefore 
the algorithmic method we have outlined in |2.4.2| may be considered reliable. 
And the use of the companion matrix (see Appendix B and C) in performing nu- 
merical calculations of the LSO parameters has been successfully demonstrated. 



2.5.2 Conversion from the BL to the spherical coordinate system. 

The foregoing analysis was performed in the BL coordinate system in which 
we suppressed the use of the, BL, subscript. To apply these estimates of the 
LSO parameters in the problem of setting the boundary conditions needed in 
modelling the evolution equations, it is necessary to convert them to the spher- 
ical coordinate system. We shall describe this conversion process, and state the 
appropriate caveats. 

Equation ([6} provides the means to convert any radial distance on an el- 
liptical orbit (that lies in the equatorial plane of the KBH) expressed in BL 
coordinates into a radial distance in spherical (or cylindrical) coordinates. But 
one cannot proceed precipitously; an elliptical orbit in the BL coordinate sys- 
tem, will be only a good approximation of an ellipse once expressed in the 
spherical coordinate system. In addition, careful consideration must be given 
to the values of <\> in their respective coordinate systems as there will be some 
important differences that will demand a more profound understanding and a 
more cautious interpretation. 

Consider the case of a test-particle in an elliptical orbit about a KBH. The 
absence of the parameter, cf>, from equation (TgI) notwithstanding; the angle, 



4>s P hericai = 4>bl T arctan j (69) 

vis. equations (|2j and Q, will force the points on the orbit that correspond 
to Rlso {spherical), Rm ax (spherical), and the position of the MBH at the focus 
of the ellipse, to be no longer collinear. Therefore the use of the values of 
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Rlso (spherical) and R M ax (spherical) to calculate e sph ericai {vis. equation ( |59| ) 
is potentially a source of error, especially for KBHs of large spin. 
We calculate: 



Rm in (spherical) — \] ^Min (BL) ^ ^ ft®) 

and 



Rm ax (spherical) — \J R\l ax (B L) + ^ ' (^) 

These two values are used to calculate the semi-major axis: 



R]\4in (spherical) ~t~ RMax (spherical) /TnA 

-^-Spherical ^ 3 \' J 

from which one can obtain 

RMin (spherical) /toA 

£ spherical ^ ■ I '"J 
A LSO (spherical) 

The latus rectum can be calculated using, 

^spherical = ^-spherical (l — £ spherical) 1 (^4) 

which is analogous to equation (p0|. The orbital frequency, W sp h er icah is ob- 



tained from equation ( 67 1 . The values of these parameters expressed in spherical 
coordinates are reported in Table [2] 

The behaviour of 4> sp hericai is not part of this study; but further investi- 
gation will be undertaken since an understanding of 4> S phericai is essential for 
properly characterising the zoom and whirl of the test-particle in its orbit. A 
diagrammatic comparison of test-particle orbits in the BL and spherical coor- 
dinate systems is shown in Figures §and[7] for a KBH of spins of S = 0.5 and 
S = 0.99 respectively. The orbit parameters are taken from Tables [i]and[4] for 
e = 0.7. One can view the shift in the value of 4> sp hericai as arising from the 
Lense-Thirring precession [33]; the orbit has a shape that can be approximated 
as an ellipse that is precessing. The orbit diagrams shown in Figures [6] and [7] 
exclude this orbital precession. 

2.5.3 LSO formulae. 

The LSO formulae we seek will be used in future work to calculate the test- 
particle orbital frequency, W, in terms of the eccentricity of the orbit, e, and 
KBH spin, S. One such relationship is already known for the SBH, i.e. 



6 + 2e (75) 
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|16l |8]; but we require additional formulae for KBH systems of various values 
of spin, and for the prograde and retrograde orbits. To this end, the algorithm 
outlined in Section |2.4.2| was used to calculate a sequence of latus rectum values, 
I, for LSOs of differing eccentricity, e, in spherical coordinates. These results 
are plotted in Figures [4] and [5] for prograde and retrograde orbits respectively, 
and each set was fit to a sixth order polynomial equation of the form, / = ct s k , 
where Ck corresponds to the coefficients to be calculated (see Tables [5] and [6} . 
The result for the SBH system is shown in each of the two figures where the 
least squares fit yielded a linear result, I = 6.00 + 2.00 e, which is consistent 
with equation (75). Such agreement is noteworthy because the least squares 
fit, based upon results previously known through the analytical and numerical 
analysis described in Section |2.4.2| corroborate the LSO relationship for the 
SBH. 

The linear approximations obtained for the KBH systems were used to cal- 
culate the LSO radii which are essential for determining the point at which an 

inspiraling CO will plunge. Although the data point pairs, (e,tj, derived for 
a particular spin became slightly nonlinear with increasing spin, the square of 
the correlation coefficient equals 1. 



3 Conclusions 

A knowledge of the relationship between the latus rectum, I, of a last stable 
orbit (LSO) and the Kerr black hole (KBH) spin, S, where S = |J| /M 2 , is es- 
sential for the calculation of the compact object (CO) orbit evolution in extreme 
black hole systems, and thus the gravitational wave energy emission. The Kerr 
metric provides the basis for the derivation of analytical relationships between 
orbital angular momentum squared (L 2 ) and the apogee of the last stable orbit 
(H-Min) for the prograde and retrograde elliptical orbits of test-particles about 
a KBH. These formulae lead directly to new and simplified representations of 
Rlso with respect to KBH spin for circular orbits, which in turn are used as a 
starting point in performing numerical analysis of elliptical LSOs. By using the 
prograde and retrograde relationships between the values of L 2 and Ruin that 
we have derived from the effective potential, an elliptical LSO can be analysed 
numerically to yield values for Rjain- The algorithm provides a foundation that 
will be generalised to include inclined orbits for which the effective potential is 
more complicated than that for the case where the orbit lies in the equatorial 
plane of the KBH. Therefore finding the relationship between Rmiu and orbital 
angular momentum becomes paramount as it allow for the methodical treatment 
of orbits of successively greater eccentricity and orbital angular momentum. 

Formulae for orbital energy, E, and the quantity, (L—SE), have led us to the 
derivation of analytical expressions for I in terms of S and orbit eccentricity. The 
LSO values obtained by using these formulae were in excellent agreement with 
those in the literature, therefore, demonstrating their validity. The usefulness 
of these analytical expressions may be found in the advantage gained in future 
theoretical and numerical investigations. These equations and the others we 
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have derived here also demonstrate the importance of using parameters that 
are normalised with respect to the KBH mass, M. 

The values of Ruin and Ruax, in Boyer-Lindquist (BL) coordinates, must be 
transformed to spherical coordinates. The I lso( spherical) an d LSO eccentricity 
(s Spherical) can then be estimated and used in the integration of the post New- 
tonian orbital evolution equations. Because we can now calculate analytically 
the esphericai and I lso (spherical) values for a range of KBH spins, (0 < S < 1; 
retrograde and prograde), it would facilitate the modelling of CO orbit evolution 
about a massive KBH. 

The companion matrix (CM) has been shown to be of great use in finding the 
roots of complicated polynomials in an analytical form. The use of the CM in 
numerical work is also encouraging, especially because one can perform various 
linear operations on the CM in order to transform the final results. 

Further investigation will be performed using the results of this work as a 
foundation. The methodologies that underlie our numerical algorithm will be 
extended to the case of inclined orbits. The radial frequency behaviour will 
also be treated by performing analytical integration of the radial path of the 
test-particle between the orbit pericentre and apocentre. The post Newtonian 
evolution equations that describe the inspiral of COs in extreme binary black 
hole systems will then be modelled. 
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Appendix A: The Inverse Kerr Metric 



The inverse Kerr metric expressed in the Boyer-Lindquist coordinates system. 
To simplify the presentation of the metric, we define the parameter: 



S = p 2 = M 2 (R 2 + S 2 cos 2 (#)). (76) 



The inverse Kerr metric is: 



^^7 



= or 1 











E(fl 2 +S 2 )+2 S 2 R-2 cos 2 (B) S 2 R 
(r 2 -2R+S 2 ) 

R 2 - 2 R + S' 2 

1 



m(r 2 -2R+S 2 ) 











-2 SR 
m(r 2 -2R+S 2 ) 




R 2 -2 R+S 2 cos 2 (0) 
(« 2 -2fl+S 2 )sin 2 (e) 



In this study, = \ , therefore, the inverse Kerr metric simplifies to the form: 



(78) 



R i +R 2 S 2 +2S 2 R 
(r 2 -2R+S 2 )r 2 






-2 



.R-2 fl+S z 
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M 2 R 2 





M(R 2 -2 R+S 2 )R 




R 2 -2R 

(fl. 2 -2fl+S 2 )M 2 fl 2 



M(fl 2 -2fl+S 2 )« 

The determinant of the Kerr metric was calculated to be, Det = — S 2 sin 2 (9). 
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Appendix B: Use of the Companion Matrix to 
Solve a Quartic Equation 

Given the task of finding the roots of a polynomial, (p(R) — 0), one might 
proceed by regarding it to be the characteristic polynomial of a matrix for 
which the eigenvalues are sought (i.e. the companion matrix) (see chapter 7 in 



p(R) = (i? 4 -12 R 3 -6 S 2 R 2 + 36 R 2 -28 S 2 R + 9 S 4 
= 0. 



(79) 



The creation of said matrix proceeds trivially to produce the companion matrix, 
M (see section 7.4.6 in |34]1: 



M = 



-9S 4 

1 28 S 2 

1 6^ 2 -36 

1 12 



(80) 



from which one may calculate the eigenvalues. These eigenvalues represent the 



solutions of equation (79 1. There are four solutions, which are (in simplified 
form) : 



where 



R 



:-! \ vZd \i Ui ; S: ' Z • 3(3 • v ! 



z 



(81) 



3 + 5 2 



We know by evaluating the solutions at S — (the Schwarzschild case) which 
two of the four solutions ought to be retained. They are: 



R 



3 + Vz± 



16S 2 



Z 



Z + 3 



(3 + 5 2 ). 



(82) 
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Appendix C: Use of the Companion Matrix to 
Find the Analytical Solution for Ilso f° r a Gen- 
eral Elliptical Orbit 

Treatment of the orbital energy, E, and the quantity, (X = L — SE), leads to 
an analytical expression for the latus rectum, I, of the last stable orbit (LSO) 
of a test-particle. An analytical form of E (see [2]) is: 



E - ^-d-,^!-^!))!-, (S3) 
In Appendix A of [2] the term X 2 has been calculated to be, 



X> = (84) 
for which the negative sign corresponds to a prograde orbit and the positive 



sign corresponds to a retrograde orbit. The functions in equation (84 1 are: 

f(7-3-£ 2N ) 2 -4,S 2 (l-£) 2 (l + £) 2 
' ' ; (85) 



; 3 

and 



l(l-3-e 2 )+S 2 (l + l + 3e 2 ) 
2^ '- ^ '-■ (86) 



/ 

and the discriminator (d = n 2 — 4/c): 



d = ) ^^l~2-2e^+S 2 (l + e) 2 ^{l(l-2 + 2e^+S 2 (l-e) 2 (^7) 
where 

c = (l-S 2 ) 2 . 

The analytical relationship between Ilso of the LSO orbit and S and e can 
be found by solving either of the following equalities: 

Xprograde (1 + e) (3 - e) = I 2 

and 
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X Retrograde (1 + £ ) (3 _ e) — ' 



(89) 



By manipulating either of the equations (|88| or (891 and removing the square 
root, one obtains the characteristic polynomial: 



P 



(0 



= [ 4 + (-4e-12)f 3 

+ (-4 S 2 e + 2 S 2 e 2 - 6 S 2 + 4 e 2 + 24 e + 36) P 

~iS 2 (1 + e) (e 2 + 7) 1 + (1 + e) 2 (-3 + e) 2 S 4 
= 0. 



(90) 



Converting the characteristic polynomial (equation ( 90 1 ) into a companion ma- 
trix yields (see section 7.4.6 in (34]): 



AI 





1 
1 
1 



-2S z e 



1J1 



-S 4 (l + e) 2 (3-e) 2 
4S 2 (l + e) (e 2 + 7) 
4e5 2 + 65 2 -4e 2 -24£ 
4(3 + e) 



36 



(91) 



The eigenvalues of equation (91 1 can be evaluated analytically and they corre- 



spond to the roots of equation (90). Two of those roots correspond to the latus 



rectum (I) of each of the prograde and retrograde test-particle orbits. One can 
also substitute the KBH spin, S, and the eccentricity of the orbit, s, into the 
companion matrix and then calculate its eigenvalues to numerically calculate 
the values of I. 

The analytical form of the eigenvalues is complicated; but a factorised form 
is presented here to illustrate how the solutions for / were identified. The four 
eigenvalues, [i — 1..4), are: 



A, ; = (3 + e) ± y/Z, 
± 



(92) 



1.6 S2 ( il £) - Z + (3 + ef + £ 2 (1 + e) (3 - e). 



By following the same reasoning as in Appendix B, we know that the solutions 
sought will each correspond to 6 + 2e when S = 0. In that case, Z a — (3 + e) 2 , 
therefore, equation (92 I simplifies to: 

A, = (3 + e)±(3 + e)±(0). (93) 



Therefore two of the eigenvalues, where Xi = 0, are excluded. 
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Figure 2: A plot of Rlso vs - L for the first and second derivatives of V with 
respect to R. 
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Figure 6: A comparison of orbits in BL and spherical coordinates for a KBH 
of spin, S — 0.5 (prograde and retrograde). The view is taken from above the 
KBH equatorial plane. Orbital precession is not included. 
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Figure 7: A comparison of orbits in BL and spherical coordinates for a KBH 
of spin, S — 0.99 (prograde and retrograde). The view is taken from above the 
KBH equatorial plane. Orbital precession is not included. 
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